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Abstract 

New approach in classification of integrable hydrodynamic chains is estabhshed. 
This is the method of the Hamiltonian hydrodynamic reductions. Simultaneously, 
this approach yields explicit Hamiltonian hydrodynamic reductions of the Hamil- 
tonian hydrodynamic chains. The concept of reducible Poisson brackets is estab- 
lished. Also this approach is useful for non-Hamiltonian hydrodynamic chains. The 
deformed Benney hydrodynamic chain is considered. 

Contents 



1 Introduction 2 

2 The Hamiltonian approach for 2-\-l quasilinear equations 4 

3 Reducible Poisson brackets and hydrodynamic chains 7 

4 The integrabihty criteria 12 

5 Classification of simplest Hamiltonian hydrodynamic chains 16 

5.1 The first integrability criterion (extra conservation law) 16 

5.2 The second integrability criterion (generating functions of conservation laws) 18 

5.3 The third integrability criterion (the Hamiltonian hydrodynamic reductions) 20 

6 Hydrodynamic chains associated with the Kupershmidt brackets 22 

6.1 The linear Hamiltonian density 24 

6.2 The Miura type and reciprocal transformations 25 

6.3 Reciprocal transformations and generating functions of conservation laws . 26 

7 Hydrodynamic chains determined by the Kupershmidt— Manin bracket 28 

7.1 The Benney hydrodynamic chain 29 

7.2 The mixed modified Benney hydrodynamic chain 29 

7.3 Reciprocal transformations and non-linear cases 31 

8 Non-Hamiltonian hydrodynamic chains 32 

9 The Zakharov hydrodynamic reductions 34 



1 



10 Conclusion 
References 



35 
36 



keywords: Poisson bracket, Hamiltonian structure, hydrodynamic chain, Liouville co- 
ordinates, reciprocal transformation, Miura transformation. Gibbons equation, Riemann 
mapping 

MSG: 35L40, 35L65, 37K10; PAGS: 02.30.J, ll.lO.E. 



1 Introduction 

This paper is devoted to classification and integrability of 2+1 quasilinear equations and 
corresponding hydrodynamic chains (see [10], [27], [28]) by the "Hamiltonian" approach 
based on the Hamiltonian formulation of integrable hydrodynamic chains and their hy- 
drodynamic reductions. Nevertheless, this Hamiltonian approach is much more universal, 
because most of explicit hydrodynamic reductions are non-Hamiltonian and this approach 
is useful as well as for non-Hamiltonian hydrodynamic chains. However, without lost of 
generality an illustration of this approach will be given on important examples connected 
with so-called M-Poisson brackets (see details in [30]). 

A hydrodynamic type system is component quasilinear system of PDE's of the first 
order (see [5]) 

ii^ = v)iu}vi, t,j = l,2,...,N. (1) 

A hydrodynamic chain is a generalization of a hydrodynamic type system on an infinite 
component case (see details below). At the same time, a hierarchy of commuting hydro- 
dynamic chains is equivalent to a family of 2+1 integrable quasilinear equations (see [7], 
[27], [28]). Hydrodynamic chains and corresponding 2+1 quasilinear equations can be 
decomposed in infinitely many ways on families of 1+1 integrable hydrodynamic type sys- 
tems (1) (see [7], [16], [28]). Nonlinear PDE's describing such hydrodynamic reductions 
have solutions parameterized by A^ arbitrary functions of a single variable (see [7], [16]). 
These nonlinear PDE's are still not solved yet (except [28]), but infinitely many particu- 
lar explicit solutions parameterized by an arbitrary number of constants been known for 
the Benney hierarchy many years ago (see [17], [18], [20], [22], [37]). Moreover, new 
hydrodynamic reductions can be found by this Hamiltonian approach (see, for instance, 
[29] and [33]). Every hydrodynamic reduction determines some particular solution of 
corresponding 2+1 quasilinear equation. If such an equation is homogeneous, then cor- 
responding solution is self-similar. Such solutions can be found by different methods (see 
[25]), based on the twistor approach, theory of the Beltrami and the Hamilton- Jacobi 
equations, the conformal mapping and the dispersionless 9— method, differential and al- 
gebraic geometry (see also [14]). An application of the generalized hodograph method 
yields most general non-self similar solutions by the nonlinear superposition principle (see 
[35]). For instance, the "symmetry constraint" method established in [3] comes from a 
theory of 2+1 integrable dispersive nonlinear PDE's such as the Kadomtsev-Petviashvili 
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(KP) or the Veselov-Novikov (VN) equations. The Hamiltonian approach is much more 
universal, because this method exists irrespective of origin of given hydrodynamic chain. 

The problem of classification of Hamiltonian chains consists of two steps. This Hamil- 
tonian approach is used to suggest some list of new Poisson brackets connected with 
hydrodynamic chains. Such classification of the Poisson brackets should be done in fur- 
ther publications (see also [30]), at least all "reducible" Poisson brackets (see definition 
below) should be described (nevertheless, we believe that any Poisson bracket is reducible 
by some appropriate "moment decomposition"). Such linear Poisson brackets have been 
introduced by I. Dorfman in [4]. The main claim of this paper is that the Hamilto- 
nian structure completely determines hierarchy of integrable hydrodynamic chains up to 
the Miura type and reciprocal transformations (This is not a true for integrable disper- 
sive and hydrodynamic type systems). Thus, a description of the Hamiltonian structures 
means a description of integrable hydrodynamic chains. In this paper without lost of 
generality we restrict our consideration on the Hamiltonian densities depending on first 
three moments A^. Obviously, they can be linear H = + f{A^,A^), quasilinear 
H = g{A^,A^)A'^ + /(A°,A^) and fully nonlinear H(A°, A^). Since the integrable 
hydrodynamic chains determined by the quasilinear and fully nonlinear Hamiltonian 
densities can be reduced (by some appropriate reciprocal transformations) to the inte- 
grable hydrodynamic chains determined by the linear Hamiltonian densities, we restrict 
our consideration on the linear case in this paper only. Two others case are just briefly 
described. 

The paper is organized in the following order. In the second section the method al- 
lowing to extract explicit Hamiltonian hydrodynamic reductions with the aid of pseudo- 
potentials is established. In the third section several reducible multi-parametric families 
of the Poisson brackets are derived. In the fourth section the integrability criteria of 
hydrodynamic chains are considered. In the fifth section the first such a family of the 
Poisson brackets is investigated by three different methods. Corresponding integrable 
hydrodynamic chains are found by the requirement of the existence of extra conserva- 
tion law creating a commuting flow (see [21], [27]), by the method of pseudopotentials 
(see [7], [36]) and separately by an application of the phenomenological Hamiltonian 
approach. All these hydrodynamic chains are connected with integrable hydrodynamic 
chains from the sixth section by a reciprocal transformation. In the sixth section hy- 
drodynamic chains associated with the Kupershmidt Poisson brackets are considered. In 
the seventh section a classification of integrable Hamiltonian hydrodynamic chains as- 
sociated with the Kupershmidt-Manin Poisson bracket is significantly simplified for the 
linear Hamiltonian density (cf. [10]). In the eighth section the method of Hamiltonian 
hydrodynamic reductions is extended on non-Hamiltonian hydrodynamic chains. Corre- 
sponding hydrodynamic type systems are non-Hamiltonian too. In the ninth section the 
modified Kupershmidt hydrodynamic chain is investigated with the aid of the Zakharov 
hydrodynamic reductions. 
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2 The Hamiltonian approach for 2+1 quasihnear equa- 
tions 



The "Hamiltonian" approach is universal and based on the "concept" of pseudo-potentials 
(see [7], [27], [36]). Let us briefly describe and illustrate this approach on the famous 
Benney hydro dynamic chain [1] 

Al = Al''' + kA'-' Al A; = 0,1,2,... (2) 

and corresponding the Khohlov-Zabolotzkaya equation 

y/\y /i ^X)X — ^ttl 

which can be obtained with the aid of two first equations from the Benney hydrodynamic 
chain and the first equation from its first commuting flow 

A^y = + + {k + l)A''Al + kA^-^Al, A; = 0, 1, 2, ... 

by eliminating of the moments A^ and A^. The method of hydrodynamic reductions devel- 
oped for these equations (see [16]) was extended for more wide class of 2+1 hydrodynamic 
type systems (see [7]), more general 2+1 quasilinear equations and hydrodynamic chains 
(see [8] and [9]) as an integrability criterion. The above equations have component 
hydrodynamic reductions written in the Riemann invariants 

ri = /i^(r)r^, = {/ {v) + A\v))tI, z = l,2,...,iV (3) 

consistent with the generating functions of conservation laws 

^^t = ^x{^ + A'), ^^y = ^.,{^ + A'^^ + A'), (4) 

where all moments A^ are some functions of the Riemann invariants r". Following [16], 
the Gibbons-Tsarev system 

can be obtained from the compatibility conditions dk{difi) = di^dkfJ'), where 

dif^ = ^-^ , (6) 
/i* — /i 

which is a consequence that /i is a conservation law density of the hydrodynamic type 
systems (3). Taking into account all moments A^ are connected with the Benney hydro- 
dynamic chain by the Laurent series 

A' A^ 

X = fi + — + — + — + ... 7 
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one can verify that the Benney hydro dynamic chain together with any its hydrodynamic 
reductions satisfy the Gibbons equation (see [15]) 

A.-/^A. = ^[/i,-a.(^ + A°)]. (8) 

We beheve that the Gibbons-Tsarev system is Darboux integrable, because plenty par- 
ticular solutions parameterized by arbitrary constants are known. However, a general 
solution still is unknown. 

In this paper the Hamiltonian approach is established. Let us choose N arbitrary 
conservation law densities /i*^'^-' = a'^, then 



,k\2 



+ A°(u) 



(9) 



The transformation a'^'(r) is invertible. Then the Gibbons-Tsarev system (in the inde- 
pendent variables a'^ on sole function only, cf. (5)) 



[a — a 



[a - a) ^ + (a - ] 



(10) 



djA^dkA^ 



la-' — a 



0, i^j^k, 



^^AO^kA'^ ' 'djA'^dkA^ ' ' ^^A'^^jAO 

can be derived from the compatibihty conditions (cf. (6)) dk{difi) = di{dkfi), where 
di = d/da'^ and 



fi — 



/i - a™ 



Assume a'^ are fiat coordinates (see [5]), then the corresponding Poisson bracket can 
be reduced (by linear transformation of field variables a'^) to the diagonal form 



(12) 



where 6'^-' is the Kronecker symbol. 

Theorem 1: The hydrodynamic type system (9) with the local Hamiltonian structure 



-d ^ 



(13) 



is a hydrodynamic reduction of the Benney hydrodynamic chain iff A^ = S^jta^' and the 
Hamiltonian density h2 = T,ek{a''Y /6 + (A^y/2. 

Proof: Let us substitute the expression A° = T^Eka^ in (11). This ODE system 



E 



-i -1 



/i — a' 



- 1 



/i — a' 

can be integrated in the implicit form 

\ = ii-^Sk ln(/i - a^), 

5 



(14) 



where A is an integration factor. Simultaneously, the equation of the Riemann surface 
satisfies the Gibbons equation (8). 

Corollary: Since = Y^Ska^, then all other moments 



A" = - 

n 



1^5^5,(a^r\ (15) 



if Y^Sk = 0. It is easy to prove by comparison (2) with (9). Then the Laurent series 
(7) reduces to so-called "waterbag" reduction (see [17]) for the equation of the Riemann 
surface (14). If 7^ the above formula is valid (see [19]), but expressions of the 
higher moment A" can be obtained by the substitution of the Laurent series (7) into the 
"deformed" above equation 

A - ^^fclnA = /i - ^efc ln(/i - a'^'), 

because the function A can be replaced by an arbitrary function A = A(A) (see (8)). 
Remark: The famous Zakharov reduction (see [37]) 

a'l = + , h^l = d^iu'^b''), ^0 = ^ 6™ (16) 

can be obtained in the same way. In comparison with the above "waterbag" case, let us 
expand the Taylor series at the vicinity of the local parameter A 

/i« = a' + Xb' + + ... 

If first 2N conservation law densities a' and fe* are independent, then the above hydrody- 
namic reduction (without the fixation A^ = S6") is given by the moment decomposition 

A'' = J^ia^b". (17) 
Suppose a* and 6* are flat coordinates, then a simplest local Hamiltonian structure is 

- d.gy:, W - d.^, (18) 

then the Hamiltonian density ha = S(a")26"72 + (A°)V2, where, indeed, A° = 

For the first time the concept of "pseudopotentials" was introduced for 2+1 quasilinear 
equation obtained by the dispersionless limit of 2+1 dispersive systems or obtained by the 
longwave (continuum) limit of 2+1 discrete systems in [36]. The meaning of the "pseu- 
dopotential" is so-called "dispersionless limit of the Lax formulation" for the integrable 
dispersive systems. This is nothing but the generating function of conservation laws for 
hydrodynamic chains (see the above example (4)). Pseudopotentials can be effectively 
found for aforementioned 2+1 hydrodynamic type systems (see [7]), for 2+1 quasilinear 
equations (see [8]) and for hydrodynamic chains (see [30]). Thus, 

1. we suppose pseudopotentials for the Hamiltonian hydrodynamic chains are given 
(cf. (4)) in the explicit form 

fi, = dj{ft; U\U',...,U''), 
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where the functions depend on a finite number of first moments A''; 

2. we look for exphcit hydrodynamic reductions in the Hamiltonian form (13), where 

N 
k=l 

The existence of the Hamiltonian density h means the existence of local Hamiltonian 
hydrodynamic reductions. The same procedure can be used for more general local Hamil- 
tonian structures (see [5]) than (13) 

dh 

al = d4g^'^], t = l,2,...,N, (19) 

where g^^ is a constant, non-degenerate and symmetric matrix, and for nonlocal Hamilto- 
nian hydrodynamic reductions (see [6]). 

Conjecture 1: Any integrable hydrodynamic chain (and associated 2+1 quasilinear 
system) admits N component hydrodynamic reductions 

= a./(a^[/i(a),f/2(a),...,[/^^(a)) 

written in the Hamiltonian form (19). 

Below we reformulate this Hamiltonian approach in the reverse direction: if the hydro- 
dynamic chain has the Hamiltonian structure, then their explicit hydrodynamic reductions 
are the Hamiltonian hydrodynamic type systems. 



3 Reducible Poisson brackets and hydrodynamic chains 

The main observation successfully utilized in this paper is that the first local Hamil- 
tonian structure (18) of the Zakharov reduction (16) of the Benney hydrodynamic chain 
(2) can be used for a constructing the first local Hamiltonian structure (see [22]) 

Al = {A\ H2} = [kA'+-~'d, + n9.A'=+"-i]^, (20) 

where (see (17)) the moments A^ = S(a"*)'^'6™. The Kupershmidt-Manin bracket (see 
(20)) 

= [A;A'=+"~i9, + n9,A'=+"^i]5(x-x'), fc, n = 0, 1, 2, ... (21) 

reduces to the canonical Poisson bracket (see (18) and details in [5], just non-zero com- 
ponent below) 

{a^ b''} = S'{x - x') (22) 

under the above moment decomposition A^ = T,{a"^)^b"^. The Benney hydrodynamic 
chain (2) is determined by the Hamiltonian H2 = | J[A'^+{A'^Y]dx, while the Kupershmidt- 
Manin bracket is associated with the momentum Hi = J A^dx and the Casimir (an- 
nihilator) Hq = / A^dx. The Zakharov reduction of the Benney hydrodynamic chain 
(16) is determined by the Hamiltonian £2 = | J[S(a'^)^6'^ -|- {Tib^Y]dx, the momentum 
hi = J Titt'^b'^dx and 2N Casimirs h(^k) = / b'^dx, h(Ar+fc) = / a'^dx. 
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In general case the local Hamiltonian structure for the hydrodynamic type system (1) 

< = {u\h}=[g^^d.^-g-Ti,ul]^^, (23) 

is determined by the Hamiltonian h = J h{u)dx and by the Dubrovin-Novikov bracket 
(a differential-geometric Poisson bracket of the first order, see [5]) 

{u%x),u^ix')} = [g'^d,-g'Xk<n^-^^'), h3 = 1,2,...,7V, (24) 

where a fiat metric 5'*-'(u) is symmetric and non-degenerate, F^^ = \g^^{dsgmk + dugms — 
dmQsk) are the Christoffel symbols. Such local Hamiltonian structure can be written via 
so-called the Liouville coordinates /l^(u), that the corresponding Poisson bracket is 

{A\x), A"(x')} = [W'="(A)9, + 9,.>V"*^(A)]5(x - x'), k,n = 1, 2, N. 

Formally, this Poisson bracket can be extended on an infinite component case, see the 
example (21). The main problem in a classification of integrable hydrodynamic chains is 
a description of such Poisson brackets 

{A\x), A^{x')} = [W'^{A)d, + 9,>V"'=(A)](5(x -x'), k,n= 1, 2, (25) 

where the coefficients W'^"(A) satisfy the Jacobi identity (see [30]) 

(26) 

Definition 1: A Poisson bracket (25) is said to be reducible if this Poisson bracket is 
equivalent to the Dubrovin-Novikov bracket (24) under the some restriction A^ = f^{u). 

Conjecture 2: Possibly, one can find appropriate constraints A'' = /'^(u) for any 
coefficients VV^'"(A), that the Poisson bracket (25) becomes reducible to the Dubrovin- 
Novikov bracket (24) in infinitely many ways. 

Example 1: The Kupershmidt-Manin bracket (21) is reducible, because at least one 
moment decomposition A^ = S(a™)'^6"* exists (see (17), (22)). Below — 1 parametric 
family of a moment decomposition is described. 

As usual we call the "hydrodynamic type system" the 1+1 quasilinear system of non- 
linear PDE's of the first order (see (1) and details in [5]). If N ^ oo, but the number 
of components v){u) is finite for every index i and each of these components is a func- 
tion of finite number of field variables u , then such hydrodynamic type systems we call 
"hydrodynamic chains" (see details in [10], [27], [28]). Also, we re-numerate equations 
of hydrodynamic chains in order that every row determined by the index i depends on 
a number of field variables larger than i. Thus, the hydrodynamic chain is an infinite 
component generalization of a hydrodynamic type system 

= 5^V;'=(A)A^, fc = l,2,..., (27) 

n=0 
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where Vj^^{A) are functions of the moments A"^, m = l,2,...Mk; and A^^ are some 
integers. 

Let us introduce the infinite set of the moments A'^ = /'^(u), where f^{u) are some 
functions. For simphcity, below we shall describe hydrodynamic type systems (1) whose 
coefficients t'j(u) depend explicitly on number N of field variables u^. 

Example 2: The Zakharov reduction (16) explicitly depends on number of the field 
variables (because A^ = T,b^), but the Benney hydrodynamic chain does not depend 
on this number. 

If the hydrodynamic type system (1) can be written via these moments A'^, whose 
coefficients Vj{A) are independent on the above number N, then such a hydrodynamic 
type system (1) we call a hydrodynamic reduction (see definition and details below) 
of the hydrodynamic chain (27). The variables are said to be "fiat" coordinates, 
the hydrodynamic type system (1) is said to be written in the canonical form (19), if the 
Christoffel symbols are vanished (see [5], [35]). 

Let us re-compute the canonical Poisson bracket (see (19)) 



via moments A^: 



{a\a^} =g'^'{x-x') (28) 

= ^^^"^^^£^^^^ - ^^^^ 

Assume the above r.h.s. can be expressed via moments A^ only, i.e. the coefficients 
>V'^"(A) in (25) are independent on number N (see the examples below). 
Definition 2: The hydrodynamic chain (27) written in the form (25) 

Al = {A^ H} = [W'^{A)d, + d,W^\A)]^, k,n = 0, 1, 2, ... (30) 

is said to be Hamiltonian. 

Remark: The above local Hamiltonian structures (30) and corresponding Poisson 
brackets still are not investigated properly. However, plenty publications are devoted to 
their particular cases (see, for instance, [4], [21], [22]), some of them will be described 
below. The classification of so-called M brackets is given in [30] . All examples presented 
in this paper belong to this class. 

Local Hamiltonian structures (30) for hydrodynamic chains we call reducible if they 
are obtained by the above procedure directly from local Hamiltonian structures (see (28)) 
written in fiat coordinates of corresponding hydrodynamic reductions. So, in this paper we 
describe a broad class of hydrodynamic chains associated with reducible Poisson brackets. 
All other Poisson brackets (30) and corresponding integrable hydrodynamic chains will 
be discussed elsewhere. 

Let us, for instance, introduce the infinite set of moments A'^ via multi-index constant 
matrices £ 

A' = 4°)a'^ + e°, A' = a' + e^^^ a' + e\ A' = ef^la'a^a' + e^^a^a' + e^^^a' + e', 
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where, obviously, all multi-index matrices e are symmetric under a permutation of their 
sub-indexes. Then 

x') = [2eWr^4^' + e^^r-ef^Wix - x'), 
x') = [44°)r"^^;;ia"«' + 4eW^-^(0)^^ 

+5«r^5«]5'(x - x') + 2[eWr"^52laV + e^^^-g-ef^a%6{x - x% ... 

The main problem in the classification of such local Hamiltonian structures for hydrody- 
namic chains is a description of all multi-index matrices e allowing to express explicitly 
the polynomials with respect to flat coordinates in right hand sides of the above Pois- 
son bracket via polynomials of moments A"'. Then such moments are analogues of the 
Liouville coordinates for an infinite dimensional space (see [5], [24], [30]), i.e. 

{A\ A'} = [{a'^^A' + + a^^)d, + d^ia'^^A' + a'^A^ + a^^)]5{x - x'), ... 

All coefficients a^^ must be independent on N. 
Example 3: Let us introduce the moments 

= 7^r^E^^(«')'''^'' ^ = 0' ±1' ±2, (31) 

where Ei, (3 and 7 are arbitrary constants. Any symmetric non-degenerate matrix g^^ can 
be written as a diagonal matrix under a linear change of the field variables a*. Thus, 
without lost of generality we assume that diagonal metric coefficients are ttj (off-diagonal 
coefficients are absent, see (12)). Then we have (see (29)) 

TV 

{A\ A^} = + (/3n + 7 - l){a''f^''+''^+^'^^~^^al]S{x - x'). 

1=1 

In general case r.h.s. can be expressed via moments A^ in infinitely many ways for 
every fixed number A^. However, in the above construction r.h.s. cannot depend on 
number A^. Since A^ is a sum of monomials of the sole degree (3k + 7, then the sum of 
monomials of another sole degree /?(/;; + n) + 2(7 — 1) must be expressed via the moments 
j^k+n+m^ where m is an integer, which is not determined yet. Then ai = l/Si and 

{A^ A"} = [{p{k + m) + l)A'^+"+™9, + (Pin + m) + 1)9,A'^+"+'"]5(x - x'), 

where 7 = l3m + 2. Since /3 7^ 0, one can re-scale all moments A^ — > jSA''. Introducing 
new parameter / = 1/(3, one can obtain the two parametric family of the Poisson brackets 

{4-,/)' 4-,/)} = [ik + m + /)Afi7+-a. + {n + m + l)d,Al;^^^-]6{x - x'), 
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{A\A'} 
{A', A'} 
{A\A^} 



dA^_^^^ dA\. 
-Tr-9 9j:^—6{x 

dA'_, dA\ 



where 



A'' = 

k + m + 2l 



N 



k-\-7n 



+2 



A;, m G Z 



Remark: If also / G Z, then the some value of running index k* 
determines the corresponding moment 



(32) 
(m + 21) 



N 



A 



-m-2l 
(m,l) 



ly^/i In a* 



i=l 



Remark: The above Poisson bracket can be obtained with the aid of another moment 
decomposition (cf. (17)) 



N 



i=l 



a) I 



where another canonical Poisson bracket (22). 

Example 4: Let us consider the nonlocal Poisson bracket (see [13]) 



{a\ a^} 



— — ea^a^ \ dr — sa^al + eald^ ^al 



x^x ^x 



6{x — x'), 



where Sij is the Kronecker symbol and e is a constant curvature. Introducing the moments 
A'' with the aid of (31), the above Poisson bracket leads to the Poisson bracket 

{A^ A^} = [{(3{k + m) + l)A'=+"+™9^ + {(3{n + m) + l)(9^A^'+"+™]5(x - x') 

-e[{(3{k + m) + 2){(3{n + m) + 2)A^'A"9^ + {(3{k + m) + 2){(3{n + m) + 1)A^A^ 



+ {P{n + m) + 2)A-A'^ - A'^d-'A:]6{x - x'), 

where 7 = /3m + 2. 

Example 5: Let us introduce the moments (cf. (31)) 

1 ^ 
pk + 7^ 

where 

f'\z) = RmU) = «0/*' + + «2/*'-' + ... + «M. 

Then we have 



N 

{A\A-} = Y,e^[f{a^)^^^-'f\a^)^4f{a^)f'^+^''f'{a^)6{x-x■). 

i=l 



Expressing r.h.s. for some special set of parameters (3 and 7 via higher moments, one can 
obtain exactly the Dorfman Poisson brackets investigated in [4]. 



11 



If, for instance (cf. (32)), 



N 



1 n—2m 



2m 



i=l 



(m — n)" 
k + 2m — n 



N 



i=l 



where p{z, g2, Qs) is the Weiershtrass elhptic function, then corresponding Poisson bracket 
(25) is determined by 

E^IA) = {4k + 10m - 6n)Af+y-^" - {k + — ^)^?2Afi;-5- - {k + m)^?3Afi;y 



Example 6: Let us introduce the moments 



N 

= / U''{a')dc 

1=1 ^ 



where 



M 



m=0 

Then we have 

AT 



(33) 



i=l 
N 



^5,[/'=+"(aO(5'(x - x') + 



n 



i=l 



k + n 



N 



i=l 



6{x — x'). 



Since (see (33)) 



^ AT M 



1=1 m=0 

then this Poisson bracket can be written in the most compact form 



N 



M 



{A',A'} = Y^eJ'[ 



x—x 



n=l 



■m=0 



These are exactly the Dorfman Poisson brackets (see [4]). 

Remark: all above Poisson brackets are obtained by the phenomenological Hamilto- 
nian approach presented in this paper. However, one can verify directly (see [30]), that, 
indeed, these Poisson brackets satisfy the Jacobi identity (26). 



4 The integrability criteria 

In this section we establish a new approach in the classification of integrable hydrodynamic 
chains and their integrability. If a hydrodynamic type system possesses the Hamiltonian 
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structure, it is not enough for integrabihty (the "integrabihty" means the existence of 
infinitely many conservation laws and commuting flows (see [35]). However, integrable 
hydrodynamic chains possessing the Hamiltonian structure can be classified. The Hamil- 
tonian hydrodynamic reductions determine a generating function of conservation laws, 
which exists for integrable hydrodynamic chains only. 

Theorem 2: If a hydrodynamic chain has a reducible Poisson bracket, then this Hamil- 
tonian hydrodynamic chain has at least one N — 1 parametric family of the Hamiltonian 
hydrodynamic reductions. 

Proof: Without lost of generality we restrict our consideration on the Kupershmidt- 
Manin bracket (21). Assume for simplicity, that the Hamiltonian is H2 = / li2{A^, A^, A^)dx. 
Then the corresponding hydrodynamic chain (20) is 

A' = (fcA^+^a. + 2d^A^+')^ + {kA'd. + d.A')^^ + kA'^-'d.^^- (34) 

This hydrodynamic chain has the moment decomposition (31), where /9 = 7 = 1. Then 
the above hydrodynamic chain reduces to the hydrodynamic type system written in the 
conservative form 

This hydrodynamic chain has at least three local conservation laws 

,0 ^ f ,adll2 .i^HsA ,1 ^ / .n5H2 ,i9H2 ,2^H2 _ 



dA^ dA^ J ' * V 9A^ dA^ 

^ f ,.,dll2,2 ,o5H2 9H2 ,1, 5H29H2 ,9H2,2x ,n^H2(9H2 

Then the hydrodynamic reduction (35) must have the same conservation laws. From the 
first of them one can obtain the following restriction Se^ = 0. Thus, N — 1 parametric 
family of hydrodynamic reductions is found. The above hydrodynamic reduction has the 
Hamiltonian form (13) while the Hamiltonian density of hydrodynamic type system \\2 
is the reduced Hamiltonian density H2 of hydrodynamic chain (34), where A^ = SemO™', 
A^ = Ee^(a'")V2, = S£:„(a'^)V3. Then one can check the identity 

JUdh2 , j / 9112 ,<9H2 , ,-,2 9H2\ , ,• 

^ da' ^ \dA^ dA^ ^ ' dA^ I 
i=i i=i ^ ^ 

= y ( ^dA' + ^dA' + ^dA'] = dU2. 
^{dA^ dA^ dA^ J ^ 

i=l ^ ' 



Thus, indeed, — 1 parametric family of hydrodynamic reductions has the Hamiltonian 
form, and are fiat coordinates. 

The main statement: if the hydrodynamic chain (34) is integrable, then above hy- 
drodynamic reductions (35) determine the generating function of conservation laws 

P^ = ^4dA^P + ^ (36) 
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by the replacement ^ p. A deformation of the Riemann surface determined by the 
equation X = X{A;p) is described by the Gibbons equation 



A, 



dA^ J dp 



Pt - d^. 



dH2 2 

dA^^ 



dU2 

OA^ 



Thus, the problem of the classification of integrable Hamiltonian hydrodynamic chains 
is reduced to the problem of the classification of generating functions of conservation laws 
(see the section "General case" in [32]). 

Here we enumerate four different tools allowing to prove an integrability of such 
hydrodynamic chains. 

1. Extra conservation law and method of pseudopotentials. For instance, suppose the 
hydrodynamic chain (34) is integrable, then the first commuting flow determined by the 
Hamiltonian density li-^lA^ , A^ , A"^ , A^) creates the corresponding generating function of 
conservation laws is 



dU 



dU 



3^2 



dA^ 



dU 



dA^ J ■ 



(37) 



The compatibility condition dt-i{pt) = dt{pt-2) leads to exactly the same system of nonlin- 
ear PDE's in involution as obtained by the tensor approach in [10]. The equations (36) 
and (37) are nothing but pseudopotentials for corresponding 2+1 quasilinear equation 
(see, for instance, (4) - the compatibility condition dt^ipt) = dtipt'i) yields the Khohlov- 
Zabolotzkaya equation, see (2) and below). Thus, we extend the concept of pseudopoten- 
tials on the theory of integrable hydrodynamic chains (see the section 5). 

2. Reciprocal transformations. The Hamiltonian density H2 can be linear H2 = A^ + 
f{A^, A^), quasilinear H2 = ^(A^, A^)A^ + f{A°, A^) and fully nonlinear H2(A°, A\ A^). 
Thus, the generating functions (36) can be separated on three sub-classes 



Pt 



Pt 



Pt 



dx 



4 



2 ^ df{A 
p + 







A' 



dA^ 



-P 



df{A',A'] 



9iA'^,A')p' 



( dg{A 
V dA^ 



dA^ 



A^) dgiA'^^A^) 
p + 



MO 



A' 



df{A^,A') ^df{A^,A'] 
-p + 



dA^ 



dU 



2 2 
P 



dH2 9H2\ 

dA^ ^ dA^J 



The difference between these sub-classes is obvious. The main claim is that the factors 
g{A^,A^) and dli2/dA'^ can be removed under some special reciprocal transformations. 
Thus, an integrability of an arbitrary {quasilinear or fully nonlinear) Hamiltonian inte- 
grable hydrodynamic chain is reduced to an integrability of a linear Hamiltonian inte- 
grable hydrodynamic chain (see the sections 6 and 7). 

3. Asymptotic at the vicinity of each singular point. Suppose we have some generating 
function of conservation laws of polynomial type (see, for instance, the above linear case) 



Pt = da: 



N 



p_ 

N 



+ aop""-' + aip^-^ + 



+ ClN-l 



14 



Then the generation function of conservation law densities p can be determined by the 
series 

P = A + Ho + ^ii + ^ + ^ + ... (38) 

If a hydrodynamic chain is integrable, then such an expansion exists (see the section 7). 

4. Hydrodynamic reductions. This is the most universal approach (see [7], [16], [32]). 
Let us consider the equation of the Riemann surface A = X{A;p) satisfying the Gibbons 
equation (cf. (8); see [32]) 

If A = const, then the Gibbons equation reduces to the generating function of conservation 
laws pt = dxip{A;p). Suppose the moments A" are functions of Riemann invariants r^; 
i.e. we seek component hydrodynamic reductions written in the diagonal form 

ri = /^'(r)r;, 

where = dip / dp\p=pi , r* = X\p=pi and distinct values p'^ are determined by the 
algebraic equation dX/dp = 0. Then a consistency of the above hydrodynamic type 
system and its generating function of conservation laws yields the extended Gibbons- 
Tsarev system (see [32]) describing all admissible functions ip{A;p) (see details below). 
This approach is adopted for the special class of the hydrodynamic reductions 

4 = (9^^(A;a^), 

which has the Hamiltonian form (19). In this case all moments can be explicitly 
expressed via flat coordinates a'^. The method of N component hydrodynamic reductions 
(see (5) in the Riemann invariants and (10) in the conservation law densities a*) leads 
to the Gibbons-Tsarev system, whose general solution is parameterized by arbitrary 
functions of a single variable. The method of A^ component Hamiltonian hydrodynamic 
reductions reduces the Gibbons— Tsarev system to the sole ODE (see for instance 
(14)). Let us consider, for instance, (10). If the hydrodynamic reduction (9) has the 
Hamiltonian structure (13), then the moment is a some function of A = Se^a". 
Thus, the system of nonlinear PDE's (10) reduces to the ODE (A^)" = 0. Thus, indeed 
A^ = Se^a" in accordance with (14). 

Thus, the Hamiltonian approach consists of following steps: 

1. A "moment decomposition" search for given Poisson bracket (see the previous 
section). 

2. A transformation to the linear Hamiltonian case. 

3. An asymptotic investigation of generating functions of conservation laws. 

4. A derivation of the extended Gibbons-Tsarev system in involution describing "in- 
tegrable" Hamiltonians and N component hydrodynamic reductions parameterized by N 
arbitrary functions of a single variable (see [32]). 
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5 Classification of simplest Hamiltonian hydrodynamic 
chains 



The Kupershmidt Poisson brackets 

{A\ A^} = [{k + ^)A''^^-^'d, + (n + 7)a.A'=+"-^]5(x - x'), (39) 

where M is an arbitrary integer, is a particular case of so-called M-brackets (see [30]) 

{v4^ A"} = + 9^.5"''=](5(x - x'), 

where the coefficients B^'"^ depend on the ffist k + n — M elements only (if /c + n < M, 
then corresponding coefficients are constants). These Poisson brackets are equivalent to 
each other under the re-numeration A^ — > A^'^^ (where s is an appropriate integer). The 
case M = 1 was considered in [22], the case M = was considered in [21], the case 
M = —1 was considered in [23]. 

In this section the classification of integrable hydrodynamic chains determined by the 
Kupershmidt Poisson brackets (see (39), M = —1 and [23], 7 7^ 0) 

Cf. = [{k + 7)C^+"+^9. + {n + 7)9.C^+"+^]^, k = OA, 2, 

where Ho= J iio{C^)dx, is given. If the Hamiltonian density Ho(C"^) is an arbitrary 
function, then the hydrodynamic chain 

Cfi = 7H[,(C°)C^^ + {k + 27)C^+^H[;(C°)C°, A; = 0, 1, 2, ... (40) 

has just the conservation law of the energy 

5.1 The first integrability criterion (extra conservation law) 

Conjecture 3: If a hydrodynamic chain has an extra conservation law, then this hydro- 
dynamic chain is integrable, i.e. this hydrodynamic chain has infinitely many conservation 
laws and commuting flows (see examples in [21]). 

Indeed, a second conservation law can be found in the form 

dAf{C')C'] = [|/'(C°)H'o(C°)(C^)2 + 7/(C0)H^(C°)C2] , 

where 

S/'H'J = rH[„ (7 + 1)/H[; = 7/'H'o. (41) 

Thus, the integrable hydrodynamic chain (see [23]) 

C^. = (^^^^(C°)^[(l - 27)C°Cri + ik + 2^)C'^'C% = 0, 1, 2, ... (42) 
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1 — 7 

is determined by the Hamiltonian density Hq = (C°) ^-^t in general case (7 7^ —1, 7 7^ 1 

1+7 

and 7 7^ 1/2)- Then the second conservation law density Hi = C^(C°) 1-27 determines 
the first commuting flow 

C% = (1 + 7)(g°)^gr + + ^7^^t^ 

1 — 27 1 — 27 

+l±X(C<0)f^[(A; + 1 + 27)C°C^+2 ^ 37(fc + 27) ^,^,^, ^ ^ 
1 — 27 1 — 27 

Remark: These hydrodynamic chains (for the integer values M = (1 — 87)/ (1 — 27)) 
were derived in [2] (see also [26]) via the Lax formulation. They are first members of 
so-called dispersionless limit of the M— dDym hierarchy. In this paper we recover the 
Hamiltonian structure for these hydrodynamic chains. We believe that the Lax formula- 
tion is connected with this Hamiltonian structure. However, this is open question. 

In particular cases: 

if 7 = 1, then (Ho = lnC°, Hi = C\C^)-^) 

C^r = ^C'r' -ik + 2)C"=-^'^, A: = 0, 1, 2, (43) 

C!^2 = 2(C°)-2C^'+2 - 2C^(C°)-3C,'+^ - 2{k + 2)C''^\C'')-''Cl 

-2(C°)-'^[(A; + 3)C°C'=+' - 3(A; + 2)C^C'=+^]C°, A: = 0, 1, 2, 
if 7 = 1/2, then (Ho = Hi = C^e^^") 

Cfi = ^e^"C,'+^ + {k + l)C'=+^e^°C°, A; = 0, 1, 2, (44) 

+36^^° [{k + 2)0'+^ + 3(A; + l)C'+^C^]Cl k = 0,1, 2, 
if 7 = — 1, then the second conservation law 

exists for any Hamiltonian density Ho(C°). Thus, in this case the hydrodynamic chain 

cf; = -h;,(c°)c^i + ik- 2)c'+'h'^{c')cI k = o,i, 2, ... 

is integrable, if the third conservation law 

atiH2(C°, C\ C^) = a,G2(C°, C\ C\ C') 
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exists. Then, 

3 



(C0)V3C2 _ l(C'0)-2/3(C'l)^ 



3C0 27(C0)2 
and the integrable hydrodynamic chain is 

Cf, = _^(cO)-4/3[3C°C,'+^ + (fc - 2)C'=+^C°], A; = 0, 1, 2, 
9 



where the Hamiltonian density is Hq = (C°) 



2/3 



5.2 The second integrability criterion (generating functions of 
conservation laws) 

Hydrodynamic reductions of the Hamiltonian chain (42) one can seek with the aid of the 
moment decomposition given in the form (31) 

' i=i 

if 7^ {l-K)/2, where K = 0,1,2,... 

Then the corresponding hydrodynamic type system 

ci. = (7-l)9,.[(c0^h[,(C°)] (46) 

has the local Hamiltonian structure (13) 

(4T) 

where the Hamiltonian ho = / ho(C°)(ia:. 
Remark: If 7 = (1 - K)/2, then 

^ ^ i=l 1=1 

For instance, if 7 = 1/2 (see (44)), then 

N N 

C = 5^e.lnc\ = - Y.e,{d)-'\ = 1,2, ... 

i=l i=l 

In such cases the compatibility of the above moment decomposition (45) with the hydro- 
dynamic chains (42) and (44) brings to the extra parametric restriction (cf. the previous 
sections) 

N 



Y^Si = 0. (4J 



1=1 
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If 7 = 1, then hydro dynamic reductions of the Hamiltonian chain (43) 
C^^ = H[,(C°)C^+^ + {k + 2)C'=+^h;;(C°)C°, A; = O, l, 2, ... 
one can seek with the aid of the moment decomposition in the degenerate form (31) 



1 ^ 

1=1 



Then the corresponding hydrodynamic type system 



01^=9.1^] (49) 



has the local Hamiltonian structure (cf. (13)) 

, _ A6ho 

where the Hamiltonian ho = / In C^dx. 

I postulate the existence of the generating function of conservation laws 

s,. = (7-l)9.[s^H[,(C°)], 7 7^1; = a. 7 = 1 (50) 

obtained from the above hydrodynamic reductions (46) and (49) by the formal replace- 
ment c* — ^ s. This is the main ansatz of this Hamiltonian approach (see (36)). 

In the same way the generating function of conservation laws for the first commuting 
flow can be found 

The substitution 2+1 quasilinear system 



in the compatibility condition (5,1)^2 = (342)^1 yields the couple of equations 

(7 + 1) a.^nN9 = 7 



2 



whose solutions are the same as solutions of the system (41). 
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5.3 The third integrabihty criterion (the Hamiltonian hydrody- 
namic reductions) 

Suppose the Hamiltonian hydro dynamic chain (40) is integrable. Then its Hamiltonian 
hydrodynamic reduction (47) must be consistent with the generating function of con- 
servation laws (50), which is a particular case of more general generating function of 
conservation laws St = dx{V{s)v) considered in details in [32]. The compatibility condi- 
tion di{dks) = dk{dis) (written via flat coordinates c*^; see (45)) yields the ODE 

7hX' = (37 - l)h[,^ 



1-7 



whose solution hg = (C°)i-27 again (see (41)), where 



diS 



/ 1 \ ^ T 



7-1 — _'i7-l \ 7 1 hg fc^) f-l — s 7-1 



27 



S7- 



Since = —di\/ds\ (see [31]), the equation of the Riemann surface A(c; s) can be found 
in quadratures 



27-1 



C — 



C^) T-1 — S 



7-1 



Taking into account the moment decomposition (45), the equation of the Riemann surface 
can be written in the form 

oo 

1-27-fc 



s 



(51) 



fc=0 



Thus, this equation of the Riemann mapping is valid for whole hydrodynamic chain (as 
well as for any its hydrodynamic reduction). Indeed, 
Theorem 3: The Gibbons equation 



7(1-7) _J_.^r..^^^ dX 



1 - 27 



9.(s— (C°)i^) 



(52) 



connects the Hamiltonian hydrodynamic chain (42) with the Riemann mapping (51) (7 ^ 
1,7 7^1/2). 

Proof: can be obtained by the direct substitution (42) and (51) in (52). 
Remark: The Gibbons equation 



e\ _dX 



CO ds 

connects the Hamiltonian hydrodynamic chain (43) with the Riemann mapping 



fcg-(fc+i)«_ 



A;=0 



The Gibbons equation 
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connects the Hamiltonian hydrodynamic chain (44) with the Riemann mapping 

oo 

A = -lns + J]CV 



,2fc 



k=0 



Let us consider two hydrodynamic chains 

k+l 



k+1 



n=0 n=0 

related by the infinitely many invertible transformations 

C'' = C''{C\C\...,C''), k = 0,1,2,... 
Theorem 4: The hydrodynamic chain 

is equivalent to the Hamiltonian hydrodynamic chain (42) 



= —L^(C')^[{1 - 27)C°C^^ + (fc + 27)C^+^C°], 

where = C^^_2^y 

Proof: The hydrodynamic chain (54) satisfies the Gibbons equation (52) 

A, - js^^{cy/-K = ^[st + (7 - mis^^icy/'^)], 

where s = p^'^ and the Riemann mapping (cf. (51)) is given by 



k=0 



2-y-k 



.k=0 



1-27 



Thus, the invertible transformations (53) can be derived from the above formula. 

Remark: The hydrodynamic chain (40) is a particular case of three parametrii 
dro dynamic chain (see [31]) 

= H[,(C°)Cf+^ + [{ak + 2)C^-+^ + pkC'' + '^kC''~^]d^H'f^{C^), 

determined by the Poisson bracket 

{C^ C^} = [Bk,nd, + d,B^,k]Six - x'). 
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where 

This Poisson bracket looks more comphcated than (39). The moment decomposition 

where the V{s) satisfies the ODE 

2VV" = aV' +(3V' + -f 

is more comphcated than (45). However, a complexity of the generating function of 
conservation laws (see [31]) is the same as in the above example (50) 

st = d4h'{A)z'is)], 

where the Hamiltonian density can be of three types h(A) = exp A, h(A) = In A, h( A) = 
A*^ and function z{s) is given in the implicit form 

dz 



6 Hydrodynamic chains associated with the Kuper- 
shmidt brackets 

The method of hydrodynamic reductions was established in [16] (for a description of 
solutions of the Benney hydrodynamic chain) and developed in [7] (see also [8]) for a 
classification of 2+1 quasilinear systems. In this paper this method is applied to classifi- 
cation of integrable hydrodynamic chains (see [21]). The second example is given by the 
Hamiltonian hydrodynamic chain 

= (i+/3)§iBr'+/3fiBi+(*+i+2^))B'« (^^'j^k+2m' (f^)^ . * = 

associated with the Kupershmidt Poisson bracket 

{5^ 5"} = [{k + f3)B'+^d, + {n + /3)9,5^-+"]5(x - x'). (55) 

If the Hamiltonian density Hi is an arbitrary, then this hydrodynamic chain has two 
conservation laws only. The conservation law of the momentum 

B% = 84(1 + 2P)B'^ + 2/35°^ - /3Hi] (56) 
and the conservation law of the energy 

2 



Hi = a 



X 
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Taking into account that this Poisson bracket is reducible by the moment decomposi- 
tion (31) (where /? ^ 7 = 2) 

= kTwY^'^ibr^'^', (57) 
^ i=i 

the Hamihonian hydrodynamic chain (see (55)) 

B^^. = [{k + /5)E^+"9. + {n + /?)a.i?'=+"]^, A; = 0, 1, 2, ... (58) 
reduces to the hydrodynamic type system 

6i.=/?a.((6^)^+"/^||j) (59) 

where Hm= Jiim{B^, B^, ■■■,B"^)dx and hm(h) = Hm|Bfc=Bfc(b)- Thus, a description of 
all admissible Hamiltonian densities can be derived from the consistency of the above 
hydrodynamic type system and the generating function of conservation laws 

q,r. = Pd^ ( q'^-/^^) . (60) 



Then the above hydrodynamic type system (59) in the Riemann invariants has the form 

In general case the Gibbons-Tsarev system can be derived from the compatibility condi- 
tions dti{dtjq) = dtj{dtiq), where di = d/dr^ (cf. (6)) 



diq 



{k + [3)[{q'flP -q^lf^] 



03*= 



and from the compatibility conditions di{dkB^) = dk{diB^), where i ^ k and n 



1,2, ...,m. 



Without lost of generality let us restrict our consideration on the Hamiltonian density 
Hi(5°,5i). Then 

_ /3 q di[dhi/dB^) + q^lPdi{dh^/dB^) 

Oiq 



13 + 1 d\ii/dB^ {q'f/P - 



The last step before a computation of the compatibility conditions di^d^q) = dk{diq) 
is the substitution the link (into the above expression for each d^q) 

_ i^ + ^Wy^'^ - {2P+l)B^S^ - 2(3B^^ ^^^ 
(/3 + 1)11- + (2/? +1)5^1^ + 2/550 ' ' 
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which is a consequence on the first equation of the hydrodynamic chain (see (56)). 

The compatibility conditions di{dkB^) = dk{diB^) imply a last restriction on the 
dependence of the Hamiltonian density hi on and B^. 

The same Gibbons-Tsarev system but written via field variables iJ' (see (59)) can be 
derived from the compatibility conditions di{djq) = dj{diq), where di = d/dV and 

. P {Y.iP^'f^^/^dkq - q^+^/^)d,{dhi/dB^) + (E h^dkq - g)g.(ghi/gE°) 

(3 + 1 [qyp -{Uy/P]d\ii/dB^ 

Before a substitution of the above expression in the compatibility conditions di{dkq) = 
dk{diq), one must compute T.(b''y~^^^'^dkq and Th^dkq- This investigation in details will 
be presented elsewhere. 

This is a particular case of more general problem: a description of all admissible 
functions ip{B^ , B^; q) (i.e. the classification of integrable hydrodynamic chains) deter- 
mined by the generating function of conservation laws qt = dxip{B^,B^]q) (see details 
in [32]). However, if the compatibility conditions di^d^q) = dk{diq) imply the extended 
Gibbons-Tsarev system describing all admissible functions ■ip{B^,B^]q) and correspond- 
ing TV component reductions parameterized by arbitrary functions of a single variable, 
if the compatibility conditions di{dkq) = dk{diq) for integrable 2+1 quasilinear systems 
like dKP (see [16] and [7]) imply the Gibbons-Tsarev system describing component 
reductions parameterized by arbitrary functions of a single variable only, the compati- 
bility conditions di{dkq) = dk{diq) for integrable hydrodynamic chains connected with the 
Kupershmidt Poisson bracket reduce the Gibbons-Tsarev system to the sole ODE only. 
It means that the advantage of this Hamiltonian hydrodynamic reductions approach in 
the application to the Hamiltonian hydrodynamic chains, whose Poisson brackets are re- 
ducible, is that the Gibbons-Tsarev system implies just a restriction on the Hamiltonian 
density. 

Integrable hydrodynamic chains associated with the Kupershmidt Poisson bracket (55) 
can be determined by the linear Hamiltonian density Hi = B^ + f{B^), the quasilinear 
Hamiltonian density Hi = g{B^)B^ + f{B^) and the fully nonlinear Hamiltonian density 
Hi(50,5i). 



6.1 The linear Hamiltonian density 

The integrability of the hydrodynamic chain determined by the linear Hamiltonian density 
U, = B^ + f{B^) 

= + j^fiB')B': + ^-^B^rmBi k = o,i, ... 

can be examined with the aid of the method of hydrodynamic reductions, where the 
generating function of conservation laws is given by (see (60)) 

= Pd, {q'+'/^ + /'(5°)g) . (61) 

In this case the Gibbons-Tsarev system can be derived from the compatibility conditions 

di{dkq) = dk{diq) and di{dkB^) = dk{diB^), where dk = d/dr^ and 



24 



The function f{B^) also can be found from the compatibihty conditions di{dkq) = dk{diq), 
where dk = d/dh^ and 



diq 



(3 + 1 f'{B^ 



The computational result is f"'{B^) = 0. Since is a momentum density (see (56)), 
one can choose f{B^) = '-/{B^y. Thus, the linear Hamiltonian density Hi = B^ + 7(5'^)^ 
determines the integrable hydrodynamic chain 



B^.=B'^^' + ^B'B'^+2^tL^B''Bl k = 0,l,..., (62) 



whose generating function of conservation laws is given by (see (61)) 

Qtr = (3d^ + 2^B\) . (63) 

This is exactly the Kupershmidt hydrodynamic chain (see [21]) up to scaling i?^' — > 
27 ^ • 

6.2 The Miura type and reciprocal transformations 

Let us consider two hydrodynamic chains 

fc+i fc+i 
= 5^F„^(C)C:, B^ = J2g'^{B)B:, k = 0,1,2,..., 

n=0 n=0 

related by the infinitely many transformations 

B'' = B\C^,C\...,C''+^), k = 0,1,2,... (64) 

Definition 3: The transformations (64) connecting B— hydrodynamic chain (right) 
and modified C— hydrodynamic chain (left) are said to be the Miura type transfor- 
mations. 

Theorem 5: The hydrodynamic chain (42) is connected with the Kupershmidt hydro- 
dynamic chain 

Bl = + (7 - 1)B'B', + {k + 2^)B'Bl k = 0,l, 2, ... (65) 
by the reciprocal transformation 

dz = Ho(CO)dx + 7H[,'(C°)CMt\ dy^ = ~J^ dt\ 

1 — 27 

where (1 - 27)5^= = c'^+i(C70)-i+(fc+i)/(i-27) (ggg (54)), 

Proof: Can be obtained by a straightforward calculation. 

Remark: The generating function of conservation laws (50) transforms in the gener- 
ating function of conservation laws for the Kupershmidt hydrodynamic chains (65) 

g, = ^9,[(l-27)g^+750g], 
7 
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where the generating function of the Miura type transformations is s = Ho(C°)q'. 

Thus, if hydrodynamic reductions of the Kupershmidt hydrodynamic chains (M— dMKP 
hierarchy, see [2] and [26]) are known, then they can be recalculated to hydrodynamic 
reductions of M— dDym hierarchy (see also details in [34]) by the above reciprocal trans- 
formation. 

The hydrodynamic chains (62), (65) are particular cases of the Kupershmidt hydro- 
dynamic chain (see [21]) 

Bl = + ^-B'B', + {k + 5)B'Bl A; = 0, 1, 2, (66) 

which are equivalent to each other under an invertible transformation for any fixed value 
of the parameter e (see the theorem below). 
Theorem 6: The Gibbons equation 



e op 



(67) 



describes a deformation of the Riemann mapping 

oo 

A = + (1 - 6)J2b''p~''-''^^\ 

where the coefficients satisfy the Kupershmidt hydrodynamic chain (66). 

6.3 Reciprocal transformations and generating functions of con- 
servation laws 

The classification of integrable hydrodynamic chains associated with the Kupershmidt 
Poisson brackets and determined by the quasilinear and fully nonlinear Hamiltonian den- 
sities will be given in separate publication. In this sub-section we restrict our consideration 
on corresponding generating functions of conservation laws. 

Theorem 7: In the quasilinear case the generating function of conservation laws 
(see (60)) 

g,i = Pd, {g{B')q'^''P + [g'{B')B^ + f'{B')]q) (68) 
can be reduced to the linear case iff 

giB^) = {B' + 5p'+'^^\ 

Proof: An arbitrary reciprocal transformation 

dz = Fdx + Gdt, dy = dt (69) 

reduces (68) to 

Py. = (3d, {g{B')F'+'/y+'/^ + [{g'{B')B' + f'{B'))F - G/P]p) , (70) 
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where the generating function of the Miura type transformations is g = Fp. Since in the 
linear case (see (63)) 

Pyi = (3d, {q'+'/^ + B'p) , (71) 

then g{B^)F^^^/^ = 1, i.e. the function F must be a some function of B^. However, F 
is a conservation law density, B^ is a conservation law density (see (56)), then F is a 
linear function with respect to B^. Thus, indeed, g^B'^) = {B'^ + 5)~(^+^/^), where S is an 
arbitrary constant. 

The first Miura type transformation is (cf. (70) and (71)) 

50 ^ _2S^±1(B' + 5)-(2+V/3)5i + (5 _ B'')f'{B'') + /(B°). 

P 

Since B^ is a conservation law density (see (63)) in the independent variables z,y, the 
r.h.s. of the expression 

B\B^ + 5) = -25^-±^(5° + 5)-(i+i/''3)5i + [{6 - B'^)f'{B^) + /(50)](fiO + 6) 

is a conservation law density in the independent variables x, t. However, this r.h.s. must 
coincide with the quasilinear Hamiltonian density Hi = g{B^)B^ + f{B^) up to the 
momentum density i?" and an arbitrary constant e. Thus, we have the ODE of the first 
order 

-25^/(5°) = [{6 - B')f'{B') + f{B')]{B' + 6), 

whose solution is 

/(5°) = (fiO - 5)2+1/^(^0 ^ 
Then the quasilinear Hamiltonian density [C, is an arbitrary constant) 

Hi = (fi° + dy^'+'/f^iB' + ^(5° - 5)2+1/^] 

creates an integrable hydrodynamic chain reducible to the Kupershmidt hydrodynamic 
chain. 

Theorem 8: In the fully nonlinear case the generating function of conservation 
laws 

,72) 

can be reduced to the linear case iff 

Hi = [Bi + (^(fi°)]WT, (73) 

where y:>(B^) is a some function. 

Proof: An arbitrary reciprocal transformation (69) reduces (72) to 



Pyi = (3d, 



^i+i//3^ 1+1//3 ^ f _ (± 
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This generating function coincides with the linear case if F^'^^^'^dHi/dB^ = 1. Since Hi 
depends on and B^, then the conservation law density F{B^,B^) must be a linear 
function with respect to the Hamiltonian density Hi up to the momentum density B^ 
and an arbitrary constant 6, which can be removed by shift z ^ z + ey, where e is an 
appropriate constant. Thus, we have the ODE 



(Hi 



i+i/g ^Hi 
dB^ 



whose solution is (73) up to an insufficient factor. The function (p{B^) can be determined 
in the same way as in the quasilinear case, but we avoid this complicated computation in 
this paper. Further details can be found in [11]. Nevertheless, we would like to emphasize 
that an application of reciprocal transformations to the quasilinear case and moreover 
to fully nonlinear case significantly simplifies computations made by the Hamiltonian 
hydrodynamic reductions method (see the beginning of this section). 



7 Hydrodynamic chains determined by the Kupershmi 
Manin bracket 

In this paper we restrict our consideration on the linear case H2 = /2 + /(A°, A^) in 
general. The hydrodynamic chain (34) 



A' = + + + 1)^' (^) + 



df 



(74) 



creates the generating function of conservation laws (36) 



P^ = ^^[Y^dA^P^dA^)^ (75) 

where p is a some function of all moments A^jA^jA^, ... and "spectral parameter" A, 
which is not introduced yet. If p — >■ 00, the corresponding fiux of the generating function 
of conservation laws (75) also tends to infinity. One can seek p as the Laurent series in 
the form (38) 

- Ho Hi H2 , , 

p = A-H„i- — (76) 

If hydrodynamic chain (74) is integrable, then the infinite series of conservation laws 



where 



(77) 

/ - 1 '^-^ - - 
dfiik = dx I Hfc+1 — ^'y^HmHfc-i- 

m=0 

H - 
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can be derived by the substitution this series (76) in (75). 

Since df/dA^ is a conservation law density, then df/dA^ must be a linear function 
of conservation law densities A^ and A^ (the Casimir density and the momentum density 
of the Kupershmidt-Manin bracket) only. Thus, 

/ = «A°Ai + /3(Ai)2 + ^(A°), 

where a and /3 are arbitrary constants (a linear term proportional A^ is removed by a 
linear change of the independent variable x — >■ x + t const), (p{A^) is not determined 
function yet (if /5 7^ 0, the constant a is unessential and can be removed by the shift 
aA^ + 2(3A^ 2/3A\ see [10] and [27]). 

Let us consider sub-cases depending on a different number of such constants. 



7.1 The Benney hydrodynamic chain 

In the simplest case H_i = const, then the generating function of conservation laws (75) 
is (the central term fi const can be removed by a linear change of the independent variable 
X X + t const) 

f^t = d.{^ + nA')). 

If the hydrodynamic chain (74) 

A'l = A^+i + kA''-^f"{A'^)Al 
is integrable, then a substitution of the Laurent series (76) 

_ Ho Hi H2 

/i-A- — 

yields the infinite series of conservation laws (77) 

fc-i 
2- 

where 



m=0 



Ho = f\A'). 

Since f\A^) is a conservation law density, then f\AP) must be a linear function of the 
conservation law density A^. Thus, one can choose f{A^) = (74^)^/2 up to unessential 
constants, which can be removed by the Galilean transformation, a shift and a scaling. 
Thus, the first classification result is precisely the Benney hydrodynamic chain (2). 



7.2 The mixed modified Benney hydrodynamic chain 

Let us consider the most general case when all constants are not vanished. In this case 
we replace the notation A'' . Taking into account that 

= a. [e^ + E'^pj , El = 3^\E' + 2E^§, + Ef^ - /]. 
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one can obtain 

Ho = 2/3^2 + 4/3^(^^)2 + iapE^E^ + y (^°)^ + 2aE^ + (1 + 2pE°)^' - 2(3^ 

from the first conservation law (77). Since the Hamiltonian density is H2 = E"^ /2 + 
aE^E^ + j3{E^y + ip{E^), then the conservation law density Hq — 4/5H2 must be a linear 
function of E^ and E^ only. It means, that Hq exists iff the function ip satisfies the ODE 
of the first order 

2 

(1 + 2PE'>)^' - + y (^°)' + jE' + S = 0, (78) 

where 7 and 6 are arbitrary constants. Integrating this equation, one can obtain the 
Hamiltonian 



^ + aE'^E' + PiE^r - ^^{Ey + (^ - iSe^Ey 



dx, 



where e is an arbitrary constant. 

Thus, the most general {mixed modified Benney) hydrodynamic chain 

2 

E^ = E^+'+{aE^+2pE')E^+{k+l)E\aE'>+2pE'),+kE''-\aE'-2ph^{Ey+{^-2pe^)E'>] 

determined by the Kupershmidt-Manin bracket and the linear Hamiltonian has the gen- 
erating function of conservation laws (instead of p we use s for this case) 

St = (y + («^° + '^(^E')s + aE' - 2f3'e\Ey + (^ - 2f3e')E'^ . 
Remark: The substitution 

fx = s + aE^ + 2(3E^ 

connects the above generating function of conservation laws with the generating function 
of conservation laws 

= 5. (y + (79) 

for the Benney hydrodynamic chain (2). 

If £ = 0, then the above hydrodynamic chain reduces to the twice modified Benney 
hydrodynamic chain 

2 

= C^+i + (aC" + 2pC^)Cl; + {k + l)C^'(ttC'° + 2pC^)^ + kC''-\aC^ + ^C°), 

Zp 

connected with the generating function of conservation laws (instead of p we use q for this 



X • 



case 



qt = d.(j + {aC' + 2PC')q + + ^C'^ 
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Remark: The substitution 

connects the above generating function of conservation laws with the generating function 
of conservation laws (79) for the Benney hydrodynamic chain (2), where 

Ho = 2(3C^ + 2aC^ + + a\C^f + Aa(3C^C^ + A(5^{C^f. 

If /? = 0, then the Hamiltonian (see (78)) 

H, = | (^ + a5°B^-^(B°)3-|(BY) dx 

determines the modified Benney hydrodynamic chain (the constant a can be removed by 
scahng aB^ ^ B^, a 0) 

2 

^ Qk+1 ^ ^^o^fc ^ ^ 1)S^S° + kB^-\aB' - y (5°)2 - 75°), 
connected with the generating function of conservation laws 

Pt = d^{j + c^B'p + aB' - ^{By - 7^°) . 
Remark: The substitution 

p = fi — aB^ 

connects the above generating function of conservation laws with the generating function 
of conservation laws 

for the Benney hydrodynamic chain (2), where = 2aB^ — 75°. 

7.3 Reciprocal transformations and non-linear cases 

In comparison with the previous section the quasilinear case H2 = g{A^, A^)A'^+f{A^, A^) 
has two sub-cases. It is easy to see by an application of the reciprocal transformation 
(69) to the generating function of conservation laws (36) 

Pt = dx{g{A^, A^)p'^ + lower order terms). 

The modified generating function of conservation laws is 

Py = dz{p^ + lower order terms), 

where 

g{A\A')F' = l. 
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If g depends on only, then the conservation law density F is a linear function of the 
Casimir density A^. Thus, g = {A^ + 5)"^. If g depends on the both moments A^ and A^, 
then the conservation law density F is a linear function of the Casimir and the momentum 
densities A^ and A^. Thus, g = {A^ + (a linear term proportional A^ is removed by 
a linear change of the independent variable x — > x + t const; see [10]). 

The fully nonlinear case is reducible to the linear case under the reciprocal transfor- 
mation (69) if 



Since F is a conservation law density depending on all three moments A'^, A^, A^, then F 
must be a linear function of the Hamiltonian density H2 up to the momentum and the 
Casimir densities, which can be removed by a shift x ^ x + et, where e is an appropriate 
constant. Thus, the Hamiltonian density 

H2= [A2 + /(A°,A^)]^/3 

determines an integrable hydrodynamic chain. The function f{A^,A^) can be found by 
the method of the Hamiltonian hydrodynamic reductions. Further details can be found 
in [12]. Nevertheless, we would like to emphasize that the "asymptotic analysis" (see the 
beginning of this section) significantly simplifies computations made by the Hamiltonian 
hydrodynamic reductions method. 



8 Non-Hamiltonian hydrodynamic chains 

The hydrodynamic chains associated with the Kupershmidt-Manin Poisson brackets (see 
(20) and (21)) 



ATT 



m+l 



possess the "waterbag moment decomposition" (15), and the corresponding hydrody- 
namic reductions are 

"m+l 



9x 



El 

n=0 



5A" 



where the Hamiltonian flm+i = J iim+i{A^,A^, ...,A^,A'^^^)dx. The hydrodynamic 
chains associated with the Kupershmidt Poisson brackets (58) also possess the waterbag 
moment decomposition (57), and the corresponding hydrodynamic reductions are given 
by (59). 

In the both above cases corresponding hydrodynamic chains linearly depend on the 
discrete variable k. 

The Hamiltonian approach presented in this paper is useful for non-Hamiltonian hy- 
drodynamic chains too. For instance, the deformed Benney hydrodynamic chain (see [10], 
cf. (2)) 

^0) 



A 



+ eA^A'^ + kA'^-'Al 



k = 0,l, 
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is not Hamiltonian, but still is integrable (at least the Hamiltonian formalism for the 
above hydrodynamic chain is unknown). Nevertheless, this hydrodynamic chain has the 
same moment decomposition (15) 



n 



+ 1 ^ 



k\n+l 



^1) 



(where Se^ = 0, cf. (48)) as the Benney hydrodynamic chain (2), but a corresponding 
generating function of conservation laws is little bit more complicated. The corresponding 
hydrodynamic type system 

(82) 



can be written in the conservative form 

1 



-(In a* - 1) a 



where 



a = exp^eu 



^3) 



Following the Hamiltonian approach, the deformed Benney hydrodynamic chain has the 
generating function of conservation laws 



Pt = 



■(Inp- 1) +eA'^)p 



^5) 



which reduces to the first equation in (4) by e — > under the substitution 

p = exp(e/i). 

Theorem 9: The deformed Benney hydrodynamic chain (80) satisfies the Gibbons 
equation (cf. (8)) 



Pt - dx 



-(Inp- 1) +eA^]p 



A,-(-lnp + eA°)A,. = |^ 
e op 

where the equation of the Riemann mapping is (7). 

Proof: can be obtained by a straightforward substitution (7) and (80) into the above 
Gibbons equation. 

Corollary: The substitution (81) in the equation of the Riemann mapping (7) yields 
again the equation of the Riemann surface (14) 



A = /i - ^£:fc ln(/i 



u 



Remark: Rational (Zakharov) reductions (of the Benney hydrodynamic chain, see 
[37]) can be obtained by a limit from the waterbag reductions (see, for instance, [3]): one 
can substitute the series u^^'^ = + v'^ /e^ + 1x1^ /{e^Y + •••5 e^'^ — > 00 to the particular case 
of the waterbag reduction 



N 



A = /i - ^^fc In 



fl — u 



k=l 
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H — u 



k • 



The same result can be obtained by a straightforward substitution the corresponding 
moment decomposition (17) 



N 



(86) 



n=l 



to the deformed Benney hydrodynamic chain. The corresponding hydrodynamic type 
system (cf. the second equation below with (82)) 



u 



k\2 



X 



can be written in the conservative form (can be obtained directly from (85) by substitution 
series p*-'^-* ^ a'' + b''\ + cf. (83) and the first equation from the system below) 



-(Ina^ - 1) a 



-\Yia' + eA^\V 



where a* = exp(£:M*) (cf. (84)) and If = ev^ exp{eu'^) 



9 The Zakharov hydrodynamic reductions 

Let us consider the integrable hydrodynamic chain (see [21]) 

= + C'>C^ + {ak + /3)C'^C° + -fkC'-'lC + ^'"^'^ {Cyi, k = 0,l, 2, ... 

This chain does not possess waterbag hydrodynamic reductions given by the moment 
decomposition (81). However, the Zakharov hydrodynamic reductions (86) exist if a 
hydrodynamic chain has a linear dependence with respect to the discrete variable k. 
Thus, the corresponding hydrodynamic reduction is 

v\ = u'vi + v'ul + C%l + pv'C°,, 
where C'^ = S(m*)'^i;*. Introducing the new field variables 

a' = (m^-7C°)^, 
the first N equations can be written in the conservative form 

7 + 1 ,„^^^2±^I±l , , ,,^oi 



V 27 + a + 1 

Thus, the generating function of conservation laws 

7+1 27+a+l ^ ^ , -1 N ^0 



Pt = d.\ . ' ^, p-^ + (7 + 1)C''P (87) 
\27 + a + 1 
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is obtained by the replacement a* p. Substituting the Taylor expansion p^''^ — > a'^ + 
Xb^ + where A is a local parameter, one obtains the second N conservation laws 



bl = d^ ((a*)^ + (7 + l)C°U* 



1-/3 



where 6* = (a*jT+it;'. 

Remark: The Kupershmidt hydrodynamic chain (66) (see [21]) has the generating 
function of conservation laws (see (67)) 



Pt = dx 



p 



SO 



e + l e 



— P 



Thus, (comparing (87) and (88)) we conclude that both above hydrodynamic chains are 
equivalent under invertible transformations = B^{C^, C^, C^), k = 0, 1, 2, where 



7 + a 



5 



/5 + 7 



7 + 1 7 + a 

Explicit expressions can be found by the comparison two series 



{p + 7C°) 



k+l ■ 



k=0^ k=0 

where the equation of the Riemann surface is given by 



a-p 

A = pT + " 



(a 



k=0 



{p + 7C0)'=+i 



Then the moments B^ can be written explicitly via the field variables a* and 6* in the 
form (cf. (17)) 



N 



B'' = (7 + a)^(a^) 7+1 b\ 



i=l 



10 Conclusion 



At this moment several very powerful approaches are appropriate in the theory of classi- 
fication of integrable hydrodynamic chains. These are the symmetry approach (see [27]), 
the method of hydrodynamic reductions (see [7], [16]), the method of pseudo-potentials 
(see [7], [36]), the tensor approach ([10]). 

However, this Hamiltonian approach (based on the concept of the Hamiltonian hy- 
drodynamic reductions) in most cases allows to check integrability and classify integrable 
hydrodynamic chains avoiding symbolic software. Moreover, this method simultaneously 
allows to construct infinitely many particular solutions parameterized by arbitrary con- 
stants in the explicit form. 
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